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DISTRIBUTED ACTUATION ALLOCATION FOR LARGE ASSEMBLIES OF 
IMPLEMENTATION UNITS 

BACKGROUND OF THE INVENTION 

1. Field of Invention 

5 [0001] This invention is directed to apparatus and methods of actuation 

allocation of a group of implementation units. 

2. Description of Related Art 

[0002] Smart matter systems are systems with a large number of elements 
such as sensors, actuators, and/or allocators. A major problem for such systems is to 

10 coordinate the actions of many elements in order to obtain a desired result or goal for 
the entire smart matter system. The desired result is achieved by controlling or 
allocating a particular action to each individual actuator to obtain a cumulative 
actuation over all the actuators. The desired individual actuation of each actuator 
can be achieved in a number of ways. However, the actuation must be accomplished 

15 in some optimal or near-optimal way while satisfying a number of constraints. 

[0003] For systems having a relatively small number of actuators, it is 
relatively easy to determine the desired individual actuation to obtain an overall 
system actuation, e.g., using an exhaustive search method. Unfortunately, the number 
of ways to produce a desired collective result increases exponentially with the number 

20 of actuators to be controlled. Thus, the problem of allocating actuation among a 

number of actuators becomes considerably more difficult as the number of actuators 
increases in the smart matter regime. Moreover, as the number of actuators increases, 
the time available to compute the allocation often decreases because of the need to 
communicate the allocation decision to the increased number of actuators. 

25 [0004] A related art method of allocating actuators involves an exhaustive 

search technique, whereby a systematic search is conducted through all possible 
actuation allocations that are able to accomplish the desired result. 

[0005] Another related art method of allocating actuation involves solving a 
mixed integer nonlinear programming problem. The allocation problem is cast as a 

30 mixed integer nonlinear programming problem over a set of continuous and integer 
variables to select the actuation for each actuator. Various methods, such as branch- 
and-bound, Bender's decomposition, and outer approximations are used to solve 



mixed integer nonlinear programming problems. These methods solve the 
optimization problem without integer constraints and then incrementally enforce the 
integer constraints. The result is that, rather than computing all of the exponential 
number of combinations that are possible, only promising branches of the search tree 
are computed with the expense of solving a nonlinear optimization problem at each 
node of the tree. 

[0006] Other related art methods to solve the constrained optimization for 
the actuator allocation use standard optimization methods such as interior point or 
active set methods. 

[0007] However, these related art methods have numerous disadvantages 
and shortfalls. 

SUMMARY OF THE INVENTION 

[0008] The related art method of allocating actuation through an exhaustive- 
search technique may operate within certain acceptable parameters for smart matter 
systems that require control of only a relatively small number of actuators, because all 
of the possible solutions can be computed in a relatively short period of time. 
However, the computational effort for exhaustive search grows exponentially with the 
number of actuators and the number of actuator states, and is therefore not feasible if 
the number of actuators increases beyond a certain limit, which is typically less than 
10 depending on the processor speed and control bandwidth desired. 

[0009] The related art method of using mixed integer nonlinear 
programming to allocate actuators also has disadvantages. The computations at each 
branch are more complex because they involve the solution of a constrained nonlinear 
programming problem at each step. Moreover, the algorithms for this method require 
full knowledge of the system in a central location, and the allocations are difficult to 
distribute among a number of processors. Thus, for methods employing nonlinear 
programming techniques, if the hardware of a system is modular in nature, extensive 
communication is required to provide the central allocator with the required 
information. 

[0010] The related art methods of allocating actuators using standard 
optimization techniques also have disadvantages. These types of methods require 
solutions that result in large deviations from optimality, because the solutions are 
restricted to binary values and numbers of actuators less than 100. These techniques 
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also require universal knowledge about the system and therefore make heavy use of 
communications, which constitutes a large cost in distributed systems. 

[001 1] As a result of the shortcomings of the above related art methods, new 
methods and apparatus for allocating actuation are needed for smart matter systems. 
5 [0012] The present invention provides methods and apparatus that 

hierarchically decompose an actuation allocation and balance the need for rapid 
computation against optimality. 

[0013] The present invention separately provides methods and apparatus that 
are able to allocate actuation for systems having differing numbers of implementation 
10 units , ranging from a few implementation units (e.g., 10) to a large number of 
implementation units (greater than 10000). 

[0014] The present invention separately provides methods and apparatus for 
allocating control of implementation units that minimize the need for communication* 
with a system controller. 
1 5 [0015] The present invention separately provides methods and apparatus that 

use continuous high level solutions to break an actuator allocation problem into 
smaller problems that can be solved in an optimal way or by breaking the smaller 
allocation problem down into yet smaller problems. 

[0016] The present invention separately provides methods and apparatus that 
20 decompose a plurality of actuators into groups of actuators and further deocompose 
the groups into smaller groups, such as modules and submodules, until the groups 
include a number of actuators that a local optimal search can be performed, wherein 
the positions of the actuators do not have to be known for the decomposition of the 
actuators and the details of a submodule are not required to be known by the module 
25 that includes the submodule. 

[0017] The present invention separately provides methods and apparatus that 
decomposes a plurality of actuators into groups of actuators wherein the grouping of 
actuators can be an arbitrary grouping. 

[0018] The present invention separately provides methods and apparatus 
30 wherein the time to compute the near optimal actuator allocation scales with the 
number of actuators. 

[0019] The present invention separately provides methods and apparatus for 
allocating actuation that are hierarchical in nature so the computation of the allocation 
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of the actuators is adaptable to the structure of the smart matter system and the smart 
matter communication system. 

[0020] The present invention separately provides methods and apparatus for 
actuator allocation that are suitable for tiled actuation. 
5 [0021] The present invention separately provides methods for actuator 

allocation that are capable of being implemented by multiprocessors. 

[0022] These and other features and advantages of this invention are 
described in or are apparent from the following detailed description of various 
exemplary embodiments of the apparatus and methods according to this invention. 
10 BRIEF DESCRIPTION OF THE DRAWINGS 

[0023] Various exemplary embodiments of systems and methods according 
to this invention will be described in detail, with reference to the following figures, 
wherein: 

Fig. 1 is a schematic of an exemplary embodiment of a many element system 
1 5 according to this invention; 

Fig. 2 is a schematic of an exemplary embodiment of an air jet paper transport 
system, according to this invention; 

Fig. 3 is a schematic of an exemplary embodiment of an array of actuators, 
according to this invention; 
20 Fig. 4 is flow chart outlining an exemplary embodiment of a method for 

allocating an overall objective for a smart matter system, according to this invention; 

Fig. 5 is a schematic of an exemplary embodiment of a physical heuristic 
grouping of an array of actuators, according to this invention; 

Fig. 6 is a schematic of an exemplary embodiment of a logical heuristic 
25 grouping of an array of actuators, according to this invention; 

Fig. 7 is a chart displaying a method of error allocation throughout a plurality 
of subdomains according to one exemplary embodiment of the present invention; 

Fig. 8 is a chart showing a pictorial representation of a sample rms combined 
error over the grid for an optimal allocation algorithm; 
30 Fig. 9 is a chart showing a grid of evaluated force/torque points (within circle); 

Fig. 10 is a chart showing the scaling of rms error and computation time for an 
optimal flat allocation algorithm; 

Fig. 1 1 is bar graph showing rms, maximum and individual errors for the 



central region of the force/torque grid for four different allocation algorithms; 

Fig. 12 is a chart showing the scaling of computation time of four different 
allocation algorithms; and 

Fig 13 is a chart showing the scaling of error of three different allocation 
5 algorithms. 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 
[0024] There are many types of many element systems that include 
implementation units, wherein an overall goal of the entire system is desired, but 
because of the aforementioned problems, it is difficult to allocate instructions to the 
10 implementation units to obtain that overall goal. The present invention provides for 
allocating that overall goal to groups of implementation units, i.e. subdomains or 
implementation modules, as well as individual implementation units. Thus, according 
to the present invention, an array of implementation units is divided into subdomains- 
where the desired goal is provided to the subdomains. Also, in accordance with the 
1 5 present invention, the subdomains can also be divided into further subdomains and the 
process can be continued until the last subdomains have a sufficiently small number 
of implementation units to perform a local optimal search. The subdomains do not 
have to align exactly with any particular implementation unit, but rather they can 
represent an arbitrary grouping of implementation units. Further, the positions or 
20 arrangements of the implementation units do not have to be known for decomposition 
of the overall goal and the details of the lower level subdomains can be hidden from 
upper level subdomains. 

[0025] Fig 1 is a schematic that shows one exemplary embodiment of a 
many element system 100 having a distributed actuation allocation, according to this 
25 invention. The system 100 can be any type of system that is used in any number of 
different types of applications, such as but not limited to an air-jet paper transport 
systems, robotic systems and the like. System 100 includes a system allocator 1 10, a 
plurality of implementation units 120, an input device 130, a controller 135 and at 
least one sensor 140. 

30 [0026] The plurality of implementation units 120 includes individual 

implementation units 150. The implementation units 150, under the control of the 
allocator 110, operate to create a desired result or action. In various exemplary 
embodiments, the plurality of implementation units 120 are an array of 
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implementation units. In various exemplary embodiments, the implementation units 
1 50 are actuators. In other various exemplary embodiments, the implementation units 
150 are air jets that act upon a piece of paper to move the paper from one position to 
another. 

5 [0027] The input device 130 communicates with the controller 135 along 

communication link 132. The input device 130 provides the controller 135 with an 
overall objective. In various exemplary embodiments the input device 130 is a 
computer or a software program configured to generate the overall objective for the 
system 100. In other various exemplary embodiments the input source 130 is, but is 

10 not limited to, a joystick, mouse, or other graphical input source read by a computer 
where the objective for the system 100 is related to the output of the graphical input 
source. In other various exemplary embodiments, the input source 130 is a higher- 
level control element or elements of a system that includes system 100. An example • 
of a higher- level control element is a job scheduler of a printing system which sends 

1 5 paper motion goals (desired trajectories) to a paper handler (an example of a system 
100). 

[00281 It should be appreciated that in various exemplary embodiments, the 
overall objective can be, but is not limited to, a trajectory of an object or in other 
various exemplary embodiments, a desired state of the implementation units. 

20 [0029] The at least one sensor 140 communicates with the controller 135 

along a communications link 142. The at least one sensor 140 monitors or detects the 
current status of the actuation of the plurality of implementation units 150 and 
provides the status information to the controller 135 along communications link 142. 
In other various exemplary embodiments, the at least one sensor 140 detects the status 

25 of an object upon which the plurality of implementation units 1 50 act. 

[0030] The controller 135 receives the overall objective from the input 
source 130. The controller 135 computes a system goal, based on the overall 
objective. In various exemplary embodiments that include at least one sensor 140, the 
controller 135 receives the current status information from the at least one sensor 140 

30 and computes instructions, based on the overall objective and the current status 

information. In various exemplary embodiments, the controller 135 is a proportional 
derivative controller. It should be appreciated that any known or later developed 
controller may be used as controller 135. 
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[0031] In various exemplary embodiments the system goal is a force 
allocation, instruction or goal for the implementation units 150 to achieve. In various 
exemplary embodiments, the system goal is a continuous value. In other various 
exemplary embodiments, the system goal is a discrete value. 
5 [0032] The allocator 110 receives the instructions from the controller 135. 

The allocator 110 implements an allocation algorithm to determine which 
implementation units should be actuated in order to achieve the overall objective for 
the system 100. The allocator 110 communicates with the plurality of implementation 
units 120 along at least one communication link 112. 

10 In various exemplary embodiments, the plurality of implementation 

units 120 are divided into groups of implementation units. The groups of 
implementation units can be further divided until the number of implementation units 
in the groups is sufficiently small enough to allow the allocator 1 10 to perform a local 
optimal allocation search. 

15 [0033] The plurality of implementation units are divided into a plurality of 

implementation unit modules 170 and 180. In various exemplary embodiments the 
implementation unit modules 170 include individual implementation units 150. In 
various exemplary embodiments the implementation unit modules 180 include at least 
one other implementation unit module 170 as a submodule. In various exemplary 

20 embodiments, more than one of the implementation unit modules 170 include the 
same individual implementation unit 150, (this arrangement is not shown in Fig. 1). 
In various exemplary embodiments, more than the implementation unit module 180 
includes the same implementation unit module 170 (this arrangement is not shown in 
Fig. 1). It is understood that the plurality of implementation units 150 may be further 

25 divided and grouped in many varying arrangements. 

[0034] In various exemplary embodiments, the allocator 110 allocates the 
overall objective to the plurality of implementation units 120 by expressing the 
overall objective in terms of solutions, module allocation, sub-goals, sub-objectives or 
sub-instructions for each of the modules 170 and 1 80. The solutions for each of the 

30 modules 170 and 180 serve as constraints for a solution for the overall objective 

allocation . The solution for each module 170 may in turn be decomposed into further 
subsolutions for submodules, thus providing hierarchical decomposition. The 
hierarchical decomposition can be continued by the allocator 110 until the 



submodules include a sufficiently small enough number of individual implementation 
units 150 so that the allocation solution for that submodule can be performed quickly 
using a conventional technique, such as but not limited to an exhaustive search 
technique. 

[0035] In various exemplary embodiments, the allocator 1 10 includes a 
plurality of allocation levels. In an exemplary embodiment, shown in Fig. 1, the 
allocator 1 10 has a top level 1 16, a middle level 117 and a bottom level 118. In 
various exemplary embodiments, the allocator 110 includes as few as one allocation 
level to as many allocation levels as necessary to communicate to all implementation 
units such that the allocation can be performed within desired time constraints. 

[0036] The top level 1 1 6 includes a module allocator 191. The middle level 

1 17 includes module allocators 192 and 193. The bottom allocation level 118 
includes module allocators 194, 195, 196 and 197. The allocation levels 1 16, 1 17 and 

1 18 are arranged in a hierarchical arrangement, whereby the module allocator 191 in 
the top level 1 16 controls the module allocators 192 and 193 of the middle level 1 17. 
The module allocator 192 of the middle allocator level 117 controls the module 
allocators 194 and 195 in the bottom allocation level. The module allocator 193 of 
the middle allocation level 117 controls the module allocators 196 and 197 of the 
bottom allocation level 1 18. The module allocators 194, 195, 196 and 197 of the 
bottom level 118 control respective groups of implementation units 170 or individual 
implementation units 150. 

[0037] It should be appreciated that the controller 135 and actuator 1 10 
shown in Fig. 1 can be implemented as portions of a suitable system, manager or 
programmed general purpose computer. Alternatively, the controller 135 and actuator 
1 10 can be implemented using an ASIC, a FPGA, a PEDL, a PLA, or a PAL, or using 
physically distinct hardware circuits, such as discrete logic elements or discrete circuit 
elements. The particular form the controller 135 and actuator 110 shown in Fig. 1 
will take is a design choice. 

[0038] The links 1 1 2, 1 32, 1 37 and 1 42 can each be any known or later 
developed device or system for connecting the array of implementation units 120, the 
goal generation device 130, and the at least one sensor 140, respectively, to the 
allocator 110, including a wire, a direct cable connection, a connection over a wide 
area network or a local area network, a connection over an intranet, a wireless 



connection, a connection over the Internet, or a connection over any other distributed 
processing network or system. In general, the links 1 12, 132, 137 and 142 can be any 
known or later developed communication system or structure usable to connect the 
respective devices. 

[0039] It should be appreciated that the plurality of implementation units 
120, the input device 130, the controller 135 and the allocator 1 10 do not have to be 
physically located together. Furthermore, it should be appreciated that the allocator 
110 and the controller 135 can be integrated into a single device. It should also be 
appreciated that the input source 130, the controller 135 and the allocator 1 10 can be 
integrated into a single device. 

[0040] Fig. 2 shows an air jet paper transport system 200, which is one 
exemplary embodiment of a many implementation unit system 100. The system 200 
includes a controller/allocator 2 1 0, which is an exemplary embodiment of a device • 
that includes a controller integrated with an allocator. In various exemplary 
embodiments the controller/allocator 210 includes hierarchical allocation levels 
similar to the allocator 1 1 0. 

[0041] The system 200 also includes an array of air jets 220 and a input 
source 230, and a plurality of sensors 240. These devices are exemplary embodiments 
of the array 120, the input device 130 and the at least one sensors 140, respectively. 

[0042] Further, the controller/allocator 210 is in communication with the 
input source 230 via communication link 232, the sensors 240 via communication link 
242 and the air jets via communication link 222. Links 222, 232 and 242 are 
exemplary embodiments of links 122, 132 and 142. 

[0043] Disposed along link 222 is at least one valve driver device 224, 
which controls the opening and closing of the individual air jet valves 250 based on 
instructions received from the allocator via link 222. The at least one valve driver 
device 224 may be any known or later developed valve driver device. 

[0044] Disposed along link 242 is at least one sensor readout 242, which 
receives signals from the sensors 240 via conventional techniques and/or devices, and 
provides the sensor signals to the controller 210 via link 242. The sensor read out 242 
may be any known or later developed device. 

[0045] An object 260 is disposed over the array of air jets 250. In this 
exemplary embodiment, the object 260 is a piece of paper. The air jets 250 will, 
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under the control of the controller/allocator 210, move the object 260 to a particular 
position or orientate the object 260 in a particular orientation. 

[0046] The sensors 240 are dispersed throughout the array 220. In various 
exemplary embodiments, the sensors 240 are CMOS image sensor bars located 
5 throughout the actuation surface. The sensors are wired to detect edge positions of the 
object 260. The spacing of the sensors 240 is such that, for typically sheet sizes of 
paper, there are always at least four points where a paper edge intersects the sensors 
240. Thus, the x, y, and angular positions of the object 260 can be determined as the 
object 200 moves over the array 220. In various exemplary embodiments, read-outs 
10 of the sensors 240 and position calculation take well less than 1 ms. 

[0047] In various exemplary embodiments, the array of air jets 220 includes 
576 angled air jets. The angled air jets are oriented in one of four different directions 
(left, right, up, and down) and distributed equally in the array 220. In this exemplary 
embodiment, the array of air jets 220 may be thought of as consisting of 36 (or 6x6) 
15 "tiles'Vactuator modules 270 of air jets. In various exemplary embodiments, each 

tile/actuator module 270 includes sixteen (or 4x4) air jets. The actuator modules 270 
can be grouped together to form larger actuation modules 280. In embodiments that 
include the large actuator modules 280, the actuator modules 270 can be thought of as 
submodule actuators of the larger module actuator 280. In accordance with the 
20 present invention, it should be appreciated that the larger actuator modules 280 could 
be grouped together to form even larger actuator modules and that the actuator 
modules 270 could be further decomposed into further actuation submodules. 

[0048] In various exemplary embodiments, individually computer- 
addressable electrostatic flap valves control the flow of air through each air jet 250. 
25 Note, the source of the air is not shown. The air from a jet impinges on the object 

260. In an exemplary embodiment, wherein the object 160 is a sheet of paper, the air 
coming from the air jets 250 makes the paper float above the array 220. In various 
exemplary embodiments, the paper floats roughly 2 mm above the exit of the air-jets 
250. Due to the viscous drag of the air against the paper, a force is applied to the 
30 paper largely confined in a small circular region around the impact point of the air on 
the paper. Forces on the order of 10" 4 Nt per jet for supply pressures of 1 kPa are 
typical. See Biegelsen, David K; Panides, Elias; Swartz, Lars E; Jackson, Warren B; 
Berlin, Andrew A. Characterization of Shear Stress Field Induced By Obliquely 
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Impinging Air Jets, MEMS-Vol. 1, Microelectromechanical Systems (MEMS) 
(ASME 1999), pp. 385-389, which is incorporated herein by reference. 

[0049] The forces of each air jet 250 add together as a linear superposition. 
Because the valves of each air jet 250 are either open or closed, each air jet 250 is 
5 essentially a binary force actuator. The entire array 220 can be viewed as a distributed 
array of binary force actuators. In various exemplary embodiments, the response time 
of an air jet 250, namely the time it takes for an air jet 250 to respond to a command, 
depends on the voltage applied to the flap valve of the air jet 250, as well as the air 
supply pressure. In various exemplary embodiments, the response time of an air jet 

10 250 is typically on the order of 1-2 ms. 

[0050] The state of object 260 can be described by six values. Particularly, 
the object 260 can have an x position, a y position, an angular position, a velocity in 
the x direction, a velocity in the y direction and an angular velocity about an axis (z • 
axis) perpendicular to the plane of the air jets. 

1 5 [0051 ] The controller/allocator 210 must choose from these multiple air j ets 

250 to apply the necessary forces to move the object 260 toward a desired overall 
objective. This arrangement illustrates that the problem of force allocation to 
actuators is an instance of the general problem of allocating a desired macro-level 
actuation to the micro-level actuators of a highly distributed system. In this 

20 exemplary embodiment, the actuators are binary air jets, which creates an integer 
programming problem. The controller/allocator 210 employs allocation methods, 
according to this invention, that solves these problems. 

[0052] As the object 260 moves across the array 220, the positions of sheet 
edges are sensed by the sensors 240. These sensor readings are used by the 

25 controller/allocator 210 to compute the x, y, and angular positions of the object 260. 
This position information over time is used to also compute the x, y and angular 
velocities of the object 260. The positions and velocities of the object 260, are the 
object state. The controller/allocator 210 computes the object state and compares it to 
a desired state given by a predetermined or dynamic trajectory, which is originally 

30 based on the overall objective from the input source 230. 

[0053] In various exemplary embodiments, the controller/allocator 210 
includes a plurality of allocation levels, each having at least one module allocator, not 
shown in Fig. 2, which are arranged in a hierarchical manner. The plurality of module 
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allocators are used to compute the desired x and y forces and a torque about the z-axis 
(perpendicular to the object 260) for each larger module actuator 280 and for each 
module actuator 270,submodule. These forces are then allocated to the individual air 
jets 250 that together apply the desired actuation forces to the object 260. In various 
5 exemplary embodiments, this control cycle is executed once every millisecond. It 
should be appreciated, that in other various exemplary embodiments, the 
controller/allocator 210 may include additional and/or separate allocation devices 
having module allocators to assist in allocating the overall objective. 

[0054] Exemplary paper transport systems are disclosed in the following: (1) 

10 Berlin, Andrew; Biegelsen, David; Cheung, Patrick; Fromherz, Markus; Goldberg, 
David; Jackson, Warren; Panides, Elias; Preas, Bryan; Reich, James; and Swartz, 
Lars, "Paper Transport Using Modulated Airjet Arrays," IS&Ts NIP 15:1999 
International Conference on Digital Printing Technologies, Orlando, Florida, October, 
1999, pp. 285-288; and (2) Biegelsen, David; Berlin, Andrew; Cheung, Patrick; 

15 Fromherz, Markus; Goldberg, David; Jackson, Warren; Preas, Bryan; Reich, James; 
and Swartz, Lars, "AirJet Paper Mover: An Example of Meso-Scale MEMS," SPIE, 
Micromachined Devices and Components VI, Santa Clara, California Vol 4176, 
September, 2000, pp. 122-129. 

[0055] Fig. 3 shows one exemplary embodiment of an array of air jets 320, 

20 which is an exemplary embodiment of the array of actuators 120. 

[0056] The array 320 includes a plurality of air jets 350. The array 320 has 
been modeled as a set of iV jets on a 2-d grid. An airjet i (i in 1 , . . . ,7V) is represented 
by its: 

position (Xrfoy); 

25 

force domains d^ and d y f, and 

forces/ x2 - and/ y ; (£; e d^,f y j e d y i). 

30 [0057] The position (*/ identified as 352, yj identified as 354) is given 

relative to the center of mass 362 of an object 360. In various exemplary 
embodiments, the center of mass 362 is approximated. However, in other various 
exemplary embodiments, the center of mass 362 is provided to the controller or the 
system is calibrated to learn the center of mass 362. In various exemplary 
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embodiments, only the actuators under the object 360 at the time of interest are 



[0058] A force domain of an individual air jet 350 may be a continuous 
interval di = [f^fmaJ or a discrete set of possible values df = {f^, f^+f^, 
. . ., f^J for some increment value f inc . In various exemplary embodiments, an air jet 
350 can only apply an x-force 390 or a y-force 395 in either a positive or a negative 
direction. In various exemplary embodiments, the air jet 350 can only be on or off 
because of a binary valve that controls the air flow through the air jet 350. Thus, for 
an x directional air-jet, the force domain is 



[0059] The equivalent holds for y-directional air-jets. In various exemplary, 
embodiments, the time for a jet valve to switch states is generally on the order of the 
control cycle, namely 1-2 milliseconds, and thus negligible relative to the motion of 
the object 360. Thus, it is assumed for the following discussion purposes that the jet 
valves can switch states (open and closed) instantaneously. 

[0060] In various exemplary embodiments, the air jets 350 may be grouped 
together to form tiles. It should be understood that each actuator module 370 may 
contain any number of individual air jets 350. The air jet module 370 may contain 
both x-directional and y-directional air jets, and thus have both an x-force component 
372 and a y-force component 374. For example, the maximum force for a module air 
jet 370 with n x-jets and n y-jets, with equal numbers of jets in positive and negative 
directions, is f MAX = n / 2 in either direction. Thus, the domains for such a module 
air-jet 370 are: 



[0061] The actuator module 370 has a position (x^ identified as 376, y x f 
identified as 378) that may be given as the centroid of the bounding box or the 
average position of air jets 350 within its bounds. 

[0062] In various exemplary embodiments, an air jet module includes a 
group of actuator modules 370, rather than individual actuators 350. An example of 
such an air jet module is identified as 380. The air jet module 380 is similar to the air 



modeled. 




djd ~ d y i — { f\lAX' - 4lAX + fma 



c...,o,... f w 



Eq. (2) 
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jet module 370 in that it can have both x and y force components, a position, and a 
force domain. 

[0063] Together, the air jets 350 in an air jet module 370 or the entire 
system or array of actuators 320 deliver an x-force F x , a y-force F y , and a z-torque T z 
to the object 360. 

[0064] In various exemplary embodiments, an actuator module * contributes 
additively with forces^/ and fa to the system's x and y-forces, respectively. 
Furthermore, the foxcQfa contributes with -yifa to the z-torque (about the center of 
mass 362 of the object 360), and the force fa contributes withx/^,/. In other words, 
the z-torque delivered by the actuator module i is t z ( = %ifa - ytfa. In various 
exemplary embodiments, the same holds for individual air jets 350, except that each 
air jet 350 delivers only either an x-force or a y-force. The total x-force F x , y-force F y , 
and z-torque T z acting on the object 360 are determined according to the following • 
linear summation model: 

r.-2£*,/„-E>,/, 

[0065] Figure 4 displays one exemplary embodiment of a method of 
hierarchically allocating a force for a many actuator system, in accordance with the 
present invention. In various exemplary embodiments, this method is employed by an 
allocator within a many actuator system. In other various exemplary embodiments, 
this method is employed by a plurality of level allocators within a many actuator 
system. 

[0066] Because the number of actuators greatly exceeds the dimensionality 
of a control command (desired goal) U, wherein U= (F x F y TJ, there are many ways 
of assigning the actuation to the actuators. In general, it is desirable to use the 
additional degrees of freedom to attain other desirable goals and thus select the most 
desirable allocation out of many possible actuator configurations. In other words, the 
force allocation problem can be viewed as a constrained optimization problem, where 
the constraints are both the domain constraints of the actuators and the relations 
between individual actuators and the total delivered forces (Eq. (2) and Eq. (3)). This 
optimization problem is solved at each control time step in order to determine which 
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actuators should be activated. Precise generation of the required action is not usually 
possible because the actuation is discrete. In general such a discrete optimization 
problem is NP-hard, i.e., actuation assignment computation scales exponentially with 
the number of actuators. On the other hand, for arrays of actuators having a large 
5 number of actuators, the actuation approaches a continuum limit, and the problem 
becomes easy to solve approximately. Therefore, a desirable property for a force 
allocation method is to solve the discrete optimization problem optimally for small 
numbers of jets, provide a near-optimal solution for intermediate numbers of jets, and 
approach the continuum for large numbers of jets, as is accomplished by the present 
10 invention. 

[0067] In various exemplary embodiments for methods of allocating forces 
according to the present invention, an optimal assignment of actuation can be obtained 
by an exhaustive search, a discrete optimization solver, or a lookup table of 
precomputed solutions for a small number of actuators (e.g., equal to or less than 10 

1 5 for the air jet example running with a loop time of 1 ms on a typical digital signal 
processor). In various exemplary embodiments for methods of allocating forces 
according to the present invention, the allocation problem is decomposed into smaller 
sub-problems in a near optimal manner using continuous solutions as approximations, 
for a large number of actuators (e.g. equal to or greater than 10 for the air jet example 

20 running with a loop time of 1 ms on a typical digital signal processor). In other 

various exemplary embodiments, each sub-problem may be further decomposed into 
yet smaller sub-problems, or, if sufficiently small, the sub-problem is solved 
optimally using conventional methods, such as, but not limited to, look-up tables. 

[0068] In various exemplary embodiments of a method for allocating a force 

25 according to the present invention, decomposition into sub-problems includes 
dividing the actuators into modules (groups of actuators) and then assigning 
responsibility to produce a required force to each module. 

[0069] The exemplary method shown in Fig. 4, begins in step S400 and 
proceeds to step S405, wherein the goal that needs to be allocated is received. In 

30 various exemplary embodiments, the goal is a force. A particular level controlled by 
either a system allocator or a module allocator is received. In various exemplary 
embodiments, the goal received is the overall goal that is received by the allocator 1 10 
from the goal generation device 130. In other various exemplary embodiments, the 
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goal received is a force allocation from an allocator 1 14 in the top level 116 received 
by an allocator 1 14 located within the middle level 117. 

[0070] In various exemplary embodiments, the top level allocation is 
identified as Allocate(i7, 7, n), wherein: U represents the overall task for the entire 
5 array of actuators and is represented as U = (F x F y T 2 ) T ; lis an identifier string or list 
of numbers, which denotes the level of the system (the top level is level 1 and an I of 
"1 .2" indicates module 2 in level 2 within the top level 1); the n denotes the number of 
modules or jets within this level (for the top level, n is 1). 

[0071] In step S410, the first module is selected, i.e. i, which is an index 
10 representing modules within a level, is set to be equal to one. For the exemplary top 
level 116, shown in Fig. 1, there is only one allocator and thus only one module. 
Accordingly, that module is selected. For the exemplary middle level 117, shown in 
Fig. 1, there are two modules, each corresponding to an allocator 114. For the 
exemplary bottom level 118, shown in Fig. 1, there are four modules, each 
1 5 corresponding to an allocator 114. 

[0072] The method proceeds to step S415, wherein the current module / is 
decomposed. To decompose a module, whether it is for the entire system, i.e. the top 
level, or for a sub level, such as the middle level, a decomposition function is used. In 
various exemplary embodiments, the decomposition function is represented as: 
20 [n^r^w,] = Decompose^, i) 

[0073] The decomposition process selects module i in level I and returns the 
number « z - of actuators or submodules, submodule position r/ = (x/ y{), and weights w/ 
for the allocation function, which are discussed further below. This decomposition 

25 function makes heuristic choices, which are discussed further below, for hierarchical 
decomposition of the module. If no (further) hierarchical decomposition is desired, 
i.e., it is not desired to split the current module into so-called sub-modules or 
subdomains, then this function decomposes module i into its nj actual jets, as would 
be done, for example, with the modules within the bottom level 118 of Fig. 1. The 

30 grouping of the plurality of actuators into actuator modules may be predetermined or 
computed dynamically based on, for example, a load-balancing criterion. 

[0074] A plurality of actuators may be grouped into actuator modules in 
varying ways. In various exemplary embodiments, the physical layout of the plurality 
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of actuators is followed, e.g., each cluster of actuators between the sensors becomes a 
module. An exemplary embodiment of a physical heuristic decomposition of an array 
of actuators is shown in Fig. 5 (i.e. "tiles"). The advantage of a physical heuristic is 
that a modular structure of the system leads to a corresponding modular structure of 
5 the allocation algorithm. In various exemplary embodiments, for an 8V 2 by 1 1 inch 
sheet of paper, which covers about 360 jets, the physical heuristic decomposes the 
area under the sheet into about 40 modules with up to about 16 actual jets in each 
module. A module entirely covered by the sheet contains 16 actual jets, with 4 jets for 
each direction. Even with discrete valves and a single-jet force of f^, such a module 
10 provides near-continuous force domains for f x i and f y f 

d* = dyi = {-4 f max , -3 f max , -2 f max , 0, 4ax> 2 f max , 3 4 4aJ 

[0075] Another exemplary grouping heuristic is a "follow the constraints" or 
a logical heuristic. An exemplary embodiment of a logical heuristic grouping of the 

15 plurality of actuators is shown in Fig. 6. For example, all x jets with the same y 

position are interchangeable with respect to the constraints in Eq. (4) which is shown 
and explained below. 

[0076] Rows of x jets with the same y position and columns of y jets with 
the same x position lead to particularly simple allocation functions within a row or 

20 column. In one exemplary embodiment, a decomposition of jets where the N modules 
are divided into iV x "x modules," each with only x-directed jets with a common y 
position yi, and N y "y modules," each with only y-directed jets with a common x 
position xj , is shown in Figure 4. With such an arrangement, the weighting features, 
( w xi w y i) = ( n d 0) f° r 3X1 x module with n x j jets and (0 n y {) for a y module containing 

25 rt y /jets. 

[0077] Grouping/decomposing at one level can greatly simplify the force 
allocation at another subordinate level. In various exemplary embodiments, for an 8V 2 
by 1 1 inch sheet of paper, the logical heuristic groups the area under the sheet into 
about 50 modules, each with up to about 12 individual jets. 
30 [0078] Referring again to Fig. 4, in step S420, the sub-goal is determined for 

the current module. In various exemplary embodiments, the sub-goal is a force. The 
sub-goal or force is represented as: 

« ( . = Solve(U, n i , r i , w i , n) 
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where U = (F x F y T Z ) T is the desired total force for all modules at the current 
level, n is the number of modules at the current level, all other inputs are those 
determined by the decomposition, and output «/ = (f K if y i tJ) T is the computed 
allocation for module i. In various exemplary embodiments, U is the overall force 
5 generated by the goal generation device 130, for example. In other various exemplary 
embodiments, U is a force allocation instruction from a module from a level that is 
higher than the level of the current module. Depending on the level in the hierarchy 
or the number of submodules in the module, step S420 will use instantiations of Eqs. 
(5) and (6), which are shown and described below. 

1 0 [0079] In step S425, it is determined whether the level that the current 

module is in is the bottom level of the entire array of actuators. The level is the 
bottom level if the module is not decomposed into further submodules. In other 
words, the current module does not allocate any forces to a module in a lower level, • 
but rather the module will provide instructions to the actuators. If the current module 

15 level is not the bottom level, then the method proceeds to step S430, wherein the sub- 
goal is allocated to modules in the next lower level. This allocation is represented as 
Allocate(« z - Li, nj). Otherwise, if in step S425 it is determined that the level that the 
current module is in is the bottom level, then the method proceeds to step S435, 
wherein the sub-goal, is assigned to individual actuators. An assignment function is 

20 represented as follows: 

u n = Assign (u ; ,I.z) 

[0080] The assignment function assigns the desired actuation uji to the 
actuators in module Li. In other words, given desired actuation ttj = (f K if y i t z i) T and a 
single actuator, the assignment function may choose/j for an x jet and^- for a y jet. 

25 For discrete domains, assignment may also apply a thresholding function to choose an 
element from the domain. For example, for an x jet with domain {0,1^}, Assign may 
return f^ if/j is greater than f max J2 and 0 otherwise. For a module with multiple 
jets, the assignment process may perform an optimal search. Since allocation is local 
to a module, it is also possible to precompute optimal assignments and reduce a force 

30 assignment to a lookup operation. 

[0081] The following methods of generating an optimal allocation and 
assignment of a desired force U given a number of modules (consisting of groups of 
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actuators or single actuators) are used for steps S430 and S435. In various exemplary 
embodiments, there are N modules with air jets in x and y directions. Each module is 
located at a position {xi,yj) {i = \,...jN) relative to the center of mass of the object, i.e., 
the module applies its force at position {xi,y(). In various exemplary embodiments, 
the jets should not work against each other while producing U (i.e., actuation should 
be minimal). Hence, a possible constrained optimization problem capturing such 
desired behavior is: 

minimize 7 Ei=i ^y- + 7 E=i ~~y 

subject to ^=Z,1./». F y=lLljy» ^ "Em*/* ( 4 > 

where f t = (f K if y i) are the allocated x and y forces for module jet i, and w s = (w x / 
w y i) are the weighting factors for each module's contribution to the x and y force. A 
large weighting factor causes the module to assume a greater role in meeting the 
constraints for forces and torque and a smaller role in minimizing the objective 
function. A small factor emphasizes the objective function. This objective function 
causes each module to minimize its actuation, while still providing the forces needed 
by the control. If activation levels are continuous, the solution to this optimization 
problem is given by: 

. F T z -F y x+F x y 
f« = — ^ — + — — ^ — ^ (y-y t ) 



fyi = W y, 



j=l V x j = i yj y j =i jy 

F v T-F v x + F x y 

-1 i 1 14 (J - x ) Eq. (5) 



^•=/ y ,(^-^)-A < (>' 1 -^) = o 

for i = 1 , . . . ,vV, where: 

i= ™yi X i - = Lm W »'^ 

2 NT* 2 

x = M * ' , j; 2 = ^ 

_L-i w * _ Em w * 

cy\=x*-x\cr)=y*-r 

[0082] The quantities x, y are the weighted average position of all the 
submodules, and <j x ,a y are the weighted standard deviations of all the x and y 



Eq. (6) 
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submodule jet positions about this average. The torque r z / is zero with respect to 

[0083] These equations provide a continuous solution to the optimization 
problem. In various exemplary embodiments wherein the domains of the jets are 
5 continuous, these equations provide an optimal solution. When implemented with 
discrete actuators, the optimal solution is approximated with a hierarchical 
decomposition as the number of jets in each module becomes large (and thus has an 
almost continuous domain). The implementation with discrete actuators will be 
discussed below. 

1 0 [0084] The overall forces required by Eqs. (5) and (6) represent a 

partitioning of the overall task U = (F x F y T z f for a large region into tasks «/ = (fa fa 

0) T for its subregions or modules. 

[0085] The solution in Eqs. (5) and (6) to the constrained optimization 

problem can be instantiated as needed for particular systems. For example, for non- 
1 5 hierarchical allocation to iV x x j ets and N y y j ets, each j et "module" includes a single 

jet and is weighted equally by (w x/ - wy) = (1 0) for an x jet and (0 1) for a y jet. The 

instantiations of Eqs. 5 and 6 are easy to derive. 

[0086] In a two-level hierarchical allocation, where the system is 

decomposed into modules consisting of individual jets, weights at the top level can be 
20 used to shift assignment of actuation depending on module sizes. For example, for the 

z'th module with x jets and n y j y jets, the weights may be selected to be (w K j w y {) = 

(« xZ - n y j). Again, the instantiations of Eqs. 5 and 6 can be derived. Allocation within 

each module is akin to non-hierarchical allocation to the jets in the module, using the 

same equations. 

25 [0087] For a module that is logically decomposed, allocation to the modules 

instantiating Eqs. (5) and (6) needs only to be computed for fa or fa for the x and y 
modules, respectively, and allocation within a module follows Eq. (7) for x modules 
and its equivalent for y modules. In various exemplary embodiments, if an array of 
air jets that has been decomposed by a logical heuristic consists only of Nx jets with 

30 the same y position, the instantiation of the allocation function for the module, 
namely, Eq. (5), would be 
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[0088] After steps S430 and S435, the method proceeds to step S440, 
wherein it is determined whether all of the modules have been evaluated, i.e. whether 
the forces for the modules have to be either further allocated or assigned to individual 
actuators. In other words, it is determined whether i is equal to n, wherein n equals 
5 the number of modules in the level. If i is not equal to n (i.e. all the modules have not 
been evaluated), the method proceeds to step S445. 

[0089] In step S445, the next module is selected. In other words, i is 
increased to The method then proceeds back through steps S415 through S440, 
until i is equal to n, the number of modules or actuators within the level. In other 

10 words, each module in the level is decomposed, the force for each module is 
determined and the force for each module is either allocated to a submodule or 
assigned to individual actuators. If there are no further modules to select, then a next 
module cannot be selected and in step S440 it will be determined that all of the 
modules have been evaluated. 

15 [0090] If in step S440, it is determined that all the modules have been 

evaluated, i.e., that i is equal to n, then the method proceeds to step S450, wherein the 
method ends. 

[0091] The foregoing exemplary embodiment has the property that the 
allocation within one module does not depend on the allocation within any other 
20 module at the same level. In other words, allocations within modules at any given 
level are completely decoupled from each other and thus can be parallelized and 
distributed among multiple processors. 

[0092] Furthermore, the distribution of computation is at the same scale as 
the actuation. This property is of obvious interest for many actuator systems with very 
25 large numbers of actuators. 

[0093] In the foregoing exemplary embodiment, each call to allocate a force 
for a module partitions the actuation to the submodules of a given module. This 
exemplary embodiment is optimal for jets with continuous domains. 

[0094] The foregoing exemplary embodiment recursively decomposes and 
30 allocates actuation to submodules within a module. In other various exemplary 

embodiments, the recursion terminates when at most « max jets remain in a module or 
the bottom level has been reached. 
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[0095] The foregoing exemplary embodiment is a recursive method. In other 
various exemplary embodiments, the method is implemented by an iterative method 
instead, as will be apparent to one skilled in the art. 

[0096] In other various exemplary embodiments, the foregoing exemplary 
allocation method further includes allocating for error in step S430 that is generated 
during the force allocation of step S430 of a previous module. The error is the 
difference between the desired goal actuation for a given module and the actually 
produced actuation. 

[0097] In one exemplary embodiment, the method of allocating the error is 
completed by a next module allocation method. The next module allocation method 
adds the force and torque errors made in allocating »,-., = (Fi A n ,Ti A n ) to the next 
module i to be allocated. In particular: 

F n = / N + (J (n l )-xF { " °) _ _ } + (1 _ w ^ AT o.-» J X; + WAF o>-» 

where: 
AF™ 

are the torque and force errors respectively made during allocation of the (z'-l)th 
subdomain and Wis an error weighting factor reflecting the relative importance of 
torque versus force errors. If W=\, then force errors are only to be considered, while 
W=0 ignores the force errors and considers only torque errors. The ith module 
attempts to correct for the error made during the earlier allocations. Similarly, when 
all the child modules/submodules have been allocated, the resulting error of the parent 
module is allocated to the next parent module who passes the error on to its children 
modules (i.e. its submodules). The order of considering the modules can be important 
with this embodiment, as well as others. For example, error allocation could be 
sequential or to a nearest neighbor. In various exemplary embodiments, an allocation 
from the rim of the object (sheet) to its center ("outside-in) is preferred. With an 
outside-in allocation, the outer most modules are allocated first. The error generated 
from this allocation is then allocated to the next outer most modules that are usually 
on the opposite side of the object. This process is continued until the center module is 
the last to be allocated. Outer modules apply the largest torques on the object. Any 
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errors made by the outer modules tend to cancel each other out, leaving smaller errors 
to be dealt with by the inner modules that can deliver smaller torques. 

[0098] In another exemplary embodiment, the error allocation method is a 
universal error allocation. In the universal error allocation method, the error 
generated is allocated to all of the remaining modules within the current level. This 
error allocation is accomplished by adjusting T( n 'V and F( n ~V by the error caused by 
allocation of the current module. The allocation of the next subdomain uses the new 
goals. Fig. 7 displays a chart that displays decomposition of a system into modules 
701, 702, and 703 for modules 1, 2, and 3 respectively, to approximate a continuous 
force represented by dot-dashed curved line 710. Each module attempts to represent 
the average continuous force density within its boundaries. For module 1, this force is 
712. The actual force delivered in this module is 706 denoted by the dashed line. The 
error therefore is given by the cross hatched square 708. The error 708 is subtracted • 
from the average continuous force 714 to give the next desired force, 718, the thick 
line in module 2. The error 708 is subtracted from 714 because the delivered force 
was too large in module 1 compared to the desired force 712. The force delivered in 
module 2 is given by 716 (dashed line), which is too small compared to the total 
desired force 718 denoted by the thick line, giving an error for module 2 given by 728. 
This error 728 is added to the desired average continuous force 722 to given the 
desired module 3 force 720 (thick line). The actual delivered force in module 3 is 724 
denoted by the dashed line. The error is therefore 726 for the combination of all three 
modules. The error scheme presented in Fig. 7 is just one of many. Another 
possibility includes assigning equal fractions of each module error to the others. 

[0099] Because of the error allocation methods, according to the present 
invention, the local allocations of force do not have to be absolutely optimal; global 
allocation is important. 

Experimental Results 
[0100] Experimental results will now be discussed to evaluate the various 
alternative force allocation algorithms in terms of their efficiency and effectiveness. 
Of particular importance for systems with large numbers of actuators is the scaling of 
their execution times and force allocation errors with the number of actuators. One 
might expect that there is a trade-off between execution time and errors for the 
different force allocation algorithms. 
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[01 01 ] For the following evaluations, a simplified system of air j ets were 
used as actuators. The system consisted of y jets that together only deliver ay-force 
and a z-torque U= [F y T Z ] T . The number of positive and negative jets Af under the 
sheet ranged from 1-100, the number of jets per module in hierarchical algorithms was 
5, and the total combined maximum force was always F max = 1 (i.e., each jet 
delivered a force of Vjv). The jets were generalized discrete jets with domains {-In, 
°> l/ N) 0-e. possible actuation states). 

[0102] There are two ways of decomposing this exemplary experimental 
system hierarchically, either by grouping jets that are next to each other into modules, 
or by defining modules as columns of multiple jets with the same x position. In either 
case, this simplified system is equivalent to controlling the angle andy-position of a 
sheet using columns of jets. 

[01 03] Each algorithm was evaluated on a two-dimensional grid of goal 
force and torque values to determine the execution time and average allocation error. 
The grid evaluation points range from -1 to 1 with step size 0.25 with a small pseudo- 
random value added to avoid symmetric cancellation. Time is the elapsed CPU time. 
All of the algorithms were implemented in Matlab 5.3 and computations for modules 
at the same level were simulated as running in parallel. The combined rms error crfor 
the desired goals U, = (F yi T zi ) and the allocated actuation = (F ya j T za i ) was 
calculated as 

where w, = w 2 = 0.5, and cr F and cr T are the standard deviations of force error 
and torque errors computed over the grid or portions of the grid. Figs. 8 and 9 show a 
task space of F y and T z grid with pseudo-random variations. Fig. 8 shows a pictorial 
representation of a sample rms combined error over the grid for an optimal algorithm. 
It shows that large force and torque values cannot be delivered together as represented 
by areas 805. Because these large, unavoidable errors dominate the rms error, the rms 
error was computed over a confined area of the grid, namely a circle 905 with range 
90% ofF max , where an algorithm should be able to allocate with a small error as 
shown in Fig. 9 (which is a grid of evaluated force/torque points (within circle)). 

[0104] In this evaluation, the four algorithms shown in Table 1 were 
compared. The four algorithms are: 1) optimal allocation using global search; 2) flat 
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(non-hierarchical) allocation using Eqs. (5) and (6) in Solve and simple thresholding 
of the result in Assign; 3) two-level hierarchical allocation using Eqs. (5) and (6) in 
Solve at both levels; and 4) two-level hierarchical allocation using Eqs. (5) and (6) at 
the top level and optimal allocation at the lower level. Optimal allocation includes a 
brute-force search through all the possibilities for the allocation that minimizes the 
rms error. 
Table 1 



Algorithms 


Optimal 


Heuristic 


Flat 


1 (global) 


2 


Hierarchical 


4 (local) 


3 



[0105] Optimal allocation provides a baseline for the average allocation 
10 error. Because the execution time of the optimal algorithm scales so poorly with the 
number of jets, the optimal algorithm was evaluated only over the range of 1 to 10 
jets. The results are shown in Fig. 10. The rms error is given by a4 //vand computation 
time scales as 0.04 10^. These scaling results are readily understood because the 
truncation error for TV jets with a force of '/preach is a uniform distribution ± X / 2 N'- 
15 When this error is averaged over the force and torque grid, the result is roughly 

YlN-Jl ' which a § rees with the observed value. While the error is small for optimal 
allocation, the exponential scaling for the computation time is unacceptable for 
systems with large numbers of actuators. 

[0106] Fig. 1 1 shows a comparison of the errors of the four algorithms in 

20 Table 1 for N = 10 jets. As the figure shows, algorithm 4 (hierarchical allocation with 
optimal allocation within the module) yields the best results among the suboptimal 
algorithms. Flat allocation is next best, while hierarchical allocation without module 
optimization has the largest errors. Comparing the hierarchical algorithm with and 
without optimal allocation indicates that most of the error occurs by truncation when 

25 assigning within the module. Thus even for a small number of jets (e.g., less than 10), 
optimal allocation yields only about a factor of two improvement in error at a cost of 1 
to 10 orders of magnitude in computation speed over the near-optimal allocation 
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schemes. It is clear that near-optimal allocation is required for the more interesting 
ranges of 10-1000 jets. 

[0107] The scaling of error and computation time as the number of jets 
ranges from 10 to 100 is now considered. While the optimal allocation cannot be 
5 computed, the scaling properties of the various near-optimal allocation algorithms can 
be compared and the scaling laws for the optimal allocation can be used. The 
computation times for the evaluated algorithms as a function of the number of jets is 
shown in Fig. 12. The hierarchical allocation algorithms exhibit a very slow increase 
with the number of jets compared with the optimal and flat algorithms, confirming the 

10 expectation that hierarchical allocations exhibit good scaling properties. The time of 
(suboptimal) flat allocation (Alg. 2) increases as N, a clear disadvantage when 
compared with hierarchical allocation. The constant penalty factor for algorithm 4 
compared with algorithm 3 reflects the time required to compute an optimal allocation 
for a small number of jets. In fact, atiV= 10, the computation time is completely 

15 dominated by the optimal allocation. This computation penalty factor could be 

removed by using precomputed look-up tables for allocation at the module level, or by 
decomposing the system such that optimal allocation within a module becomes trivial. 
In the evaluated system, the latter can be achieved by aggregating modules from 
columns of y jets with the same x position, as indicated above. 

20 [0108] The rms error for the evaluated algorithms as a function of the 

number of jets is shown in Fig. 13. As is evident, the error does not decrease 
significantly as the number of jets increases, even with optimal module allocation 
(algorithm 4). This behavior is a consequence of the fact that errors from individual 
module and jet allocations are not compensated by other allocations. More 

25 particularly, with a jet force of V^y and optimal module allocation, each module makes 
an error of up to With 5 jets per module, the number of modules is M= N I 5 , and 
thus the total force error can be as large as M X I 7 N= N / 5 7^= 0.1 independent of the 
number of jets. Because the maximum error tends to dominate the rms average, all the 
algorithms cluster near this value. Even with optimal module allocation, only the jets 

30 within a module are able to compensate for each other's error. However, from the 

maximum error formula M V^, It is apparent that if the number of modules M is kept 
constant as the number of jets Nis increased (i.e., the number of jets per module is 
increased), the error scales roughly as M x 1 /^. 
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[0109] The experimental results presented above can be summarized as 
follows. There is a trade-off between error and computation time for the allocation 
problem for large numbers of actuators. Globally optimal allocation (algorithm 1), 
while producing the best errors, scales poorly with increasing numbers of jets. In fact, 
unless optimal allocation can be precomputed or computed easily (as in the 
row/column decomposition), the computational time scales as 2^ for TV jets. On the 
other hand, flat allocation (algorithm 2), or simply breaking the problem into modules 
that are not allocated optimally (algorithm 3), results in relatively large allocation 
errors. In fact, as shown above, the dominating rms error is on the order of one half 
the maximum total force. 

[0110] The best compromise is a hierarchical decomposition into modules 
with optimal allocation within a module (algorithm 4). Given a scaling of the error as 
as shown above, the larger the module for which the optimal allocation can be • 
computed in real-time, or precomputed and stored, and thus the smaller the number of 
modules M, the better the result. At the same time, given N/ M j e ts per module, the 
optimal allocation time per module scales with 2^ , which suggests a trade-off 
between module size and number of modules for a particular TV in order to decrease 
error while keeping computation time within acceptable bounds. The relation between 
error and computation time for all the algorithms, optimal, flat, and hierarchical, can 
be summarized with the relation 

log(r)x£ccl (8) 

where Tis the computation time and E is the fractional error, the error as a fraction of 
the full range. This relationship applies for problems where, given a local minimum, a 
lower minimum can be found almost anywhere and is only weak correlated with the 
previous minimum, and/or where the algorithm does little better than chance to go 
from one local minimum to finding the next. This appears to be the case for the 
combined optimization of force and torque from discrete jets. 

[0111] However, if the algorithm can find better near-optimal solutions 
determined in constant time (i.e., precomputed or computed trivially as in the 
row/column decomposition), and if module allocation can be distributed, then even 
superlinear scaling can be achieved, as computation time is dominated by the 
allocation to M modules at the top level. 
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[0112] These results are summarized in Table 2. Note that, for hierarchical 
algorithms, the computation time is the sum of computation times at the top and 
module levels. Assuming P processors available for module allocation, the latter can 
further be divided by P (consider in particular the case where P = M, i.e., there is one 
processor per module). For reference, the communication load is included, which 
scales with the number of calls between the levels and to the actuators. 
Table 2 



Algorithms 


1 


2 


3 


4 


4 (cto) 


Rms error 


1 


1 


1 


M 


M 




IN 


2 


2 


IN 


2N 


Computation time 


2 N 


N 


M + — 
P 


N 

M + 

P 


M 

M + — 
P 


Communication load 


N 


N 


M + N 


M + N 


M + N 



[0113] Table 2 shows scaling laws of error, computation time, and 
communication time for the algorithms that were evaluated, as well as for algorithm 4 
with a constant-time module optimization (cto). The communication load is the 
number of messages passed between computational nodes and from computational 
nodes to the actuators. 

[01 14] The hybrid hierarchical-optimal algorithms presented work well for 
all ranges of numbers of actuators. For small numbers of actuators, the upper-level 
continuous allocation can be allocated trivially to a small number or even one module. 
The jets within a module are assigned optimally. For very large numbers of actuators 
(e.g., > 1000), the number of jets within a module can be large, so the continuous 
allocation (which is optimal in the continuum limit) works very well. The lower level 
optimal allocation can be done either recursively to lower levels of modules or - 
because of the large numbers - through crude but rapid assignment. The algorithm 
works well in the medium range of 10<N<1000 as well. The asymptotic optimality of 
the continuum solution is used to guide the optimization to near-optimal smaller 
problems which are then solved using optimum methods. Thus, unlike other 
algorithms, this hybrid hiearchical optimal allocation works well in all regimes, and 
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the system does not need to be redesigned as the number of actuators varies for 
different applications. 

[0115] The hybrid hierarchical allocation of the present invention has been 
applied to an actual micro-electrical-mechanical system (MEMS) and found to work 
well. An exemplary embodiment of the present invention performs the allocation in < 
0.4 ms for row/column allocation on a single Analog Devices SHARC 40MHz DSP 
using algorithm 4 with constant time optimal allocation at the lower level. An 8'/ 2 by 
1 1 inch sheet in an exemplary embodiment of the present invention covers about N = 
360 jets, and row/column decomposition results in about M= 40 to 50 modules with 
an average of 7 to 9 jets per module and a maximum of about 1 2 jets per module 
(depending on sheet position). The available jet force is about 0.1 mNt per jet. The 
average error is M / 2 f max = 20 to 25 times f max , and the time is 40 to 50 times a 
computation for a module allocation, which is dominated by maintaining the data 
structures. The resulting position error has been 30 um rms and 10" 5 radian rms angle 
error. 

[01 16] This is a tradeoff between the error and the computation time for 
methods of actuator allocation. Basically the log of the computation times the error is 
constant for hierarchical, optimal allocation, and hybrid-optimal hierarchical 
algorithms. Thus, a small decrease in error requires a large increase in computation 
time. However, because the allocation calculation is performed repetitively, this 
computational effort can be precomputed once. Thus, the hybrid hierarchical optimal 
allocation scheme had the best tradeoff between error and computation time. 

[0117] The hybrid hierarchical-optimal algorithms of the present invention 
work well for all numbers of actuators from small traditional systems to many 
actuator systems conceived with thousands of actuators. These hybrid algorithms 
represent a general approach to high-dimensional problems. In other words, 
exemplary embodiments of the present invention use a continuous solution to break 
the problem up into a number of smaller self-similar problems that can either be 
further reduced in size or solved using optimal solutions. This present invention also 
has general application in other many degree of freedom systems such as robotic 
systems, such as, but not limited to, temperature-sensitive micro-robots (i.e. "smart 
dust"), and large-scale factories. 
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[01 18] While this invention has been described in conjunction with the 
specific embodiments outlined above, it is evident that many alternatives, 
modifications and variations will be apparent to those skilled in the art. Accordingly, 
the preferred embodiments of the invention, as set forth above, are intended to be 
illustrative, not limiting. Various changes may be made without departing from the 
spirit and scope of this invention. 



